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Abstract 

A new construction is given for obtaining i?-matrices which solve the McGuire- 
Yang-Baxter equation in such a way that the spectral parameters do not possess the 
difference property. A discussion of the derivation of the supersymmetric U model 
for correlated electrons is given in this context, such that applied chemical potential 
and magentic field terms can be coupled arbitrarily. As a limiting case the Bariev 
model is obtained. 
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1. Introduction 



The study of integrable models which are constructed through the Quantum Inverse 
Scattering Method (QISM) makes available many significant and non-perturbative results 
as these models are generally solvable through Bethe ansatz methods. The quantum 
algebras |], || and their natural Z 2 -graded analogues quantum superalgebras || |J provide 
us with the principle examples of quasi-triangular Hopf (super) algebras which allow us 
to systematically construct solutions of the McGuire- Yang-Baxter (MYB) equation which 
are central to the QISM approach. The solutions of the MYB equation are referred to as 
i?-matrices. The approach of the QISM is to first build a family of commuting transfer 
matrices from the i?-matrices and then define the Hamiltonian operator as the logarithmic 
derivative of the transfer matrix evaluated at the shift point. The shift point is that value 
of the spectral parameter such that the .R-matrix behaves as the permutation operator. 
Those i?-matrices that do have a shift point are called regular and it is this property which 
ensures that this definition gives a global Hamiltonian which is a sum of local two-site 
Hamiltonians. 

The structure of the affine quantum (super) algebras is such that solutions of the MYB 
equation obtained automatically possess the difference property in the spectral parameter 
as a result of solving Jimbo's equations || |J . On the otherhand, there exist other solutions 
of the MYB equation which do not have the difference property such as Shastry's solution 
for the Hubbard model 0. Another model which has been derived through the QISM 
is the Bariev model || which was introduced as a solvable generalization of Hirsch's 
hole superconductivity model ||. In independent works by Zhou [|I(J and Shiroishi and 
Wadati [|ll[] integrability of the Bariev model was established by obtaining an appropriate 
solution of the MYB equation. Though the solutions of |TD[ and [|TTJ are different, they 
both share the feature that they do not have the difference property. Energy spectra 
obtained through the algebraic Bethe ansatz have been studied in [TjJ . 

Yet another model which is of interest in the study of correlated electron systems is 
the super symmetric U model |14| and importantly in the context of the work discussed 
here its anisotropic generalization |HJ which can be derived in the framework of the 
QISM using an i?-matrix associated with one-parameter family of minimal typical repre- 
sentations of the quantum superalgebra U q (gl (2\1). As such, the R- matrix used to derive 
this model does have the difference property. In a particular limit, the anisotropic super- 
symmetric U model reduces to the Bariev model with the addition of a divergent chemical 
potential term. The fact that these two models are related yet have been individually 
derived from i?-matrices of differing character with respect to their spectral parameter 
dependence has been somewhat mysterious. 

In this paper a construction which relates models which differ only by field terms 
such as the chemical potential will be given and a connection between i?-matrices of the 
difference property type and those without the difference property will be established. 
The connection lies in the use of twisting of the algebraic i?-matrices by a suitable twistor 
which retains the quasi-triangular Hopf-algebra. The twisting construction was developed 
T7| and gave rise to the more general notion of quasi-triangular quasi-Hopf algebras. 
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However, provided the twistor satisfies appropriate properties known as the twisted 2- 
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cocycle condition, the twisted structure can also be of the Hopf algebra type. Such twisting 
operations have lead to some significant developments such as the construction of elliptic 
solutions of the MYB equation |T8| , |19| and the construction of multiparameteric integrable 
. The latter were based on a particular type of twistor due to Reshetikhin 
a generalized form will be used which was gven by Engeldinger and Kempf 



systems |20 
ill. Here, 



. By a certain parametrization, the regularity property of the i?-matrix will be shown 



to still hold which allows the construction of an integrable one-dimensional model with 
two-site interactions. This type of parametrization is adopted from the method used in 
23| to obtain an extended region of integrability for the supersymmetric U model. 

As an example, the aniostropic supersymmetric U model will be constructed with 
arbitrary chemical potential and magentic field term. In an appropriate limit the Bariev 
model is recovered. However, an attempt to obtain an i?-matrix solution for the Bariev 
model proves fruitless and still leaves open the question of the origins of the solutions 
obtained in [jlO, 11 in terms of an underlying algebraic structure. 

Finally, it worthwhile to observe that the construction employed here may be under- 
stood in terms of representations of coloured quantum algebras |24| . (These are not to be 
confused with colour algebras in the sense of P5| .) In fact, the relation between coloured 
quantum algebras and twisting procedures has been discussed by Chakrabarti and Jagan- 
nathan f26fl , however it appears that their role in relation to integrable systems has not 
yet been addressed in the literature. 



2. Quantum Inverse Scattering Method 

Let R(u), R(u, v) G End V <8> V give a solution of the MYB equation 

R 12 (u, v)R 13 (u)R 23 (v) = R 23 (v)R 13 (u)R 12 (u, v). (1) 

For full generality, we consider the cases when V denotes a Z 2 -graded vector space. In such 
instances it is necessary to impose the following rule for the tensor product multiplication 
of matrices: 

(a ®b)(c®d) = (-l)^ac <8> bd (2) 

for matrices a, b, c, d of homogeneous degree. The symbol (a) e Z 2 denotes the degree of 
the matrix a. The transfer matrix is defined 

t(u) = str (Ron(u)R ( N - 1 )(u)....R 01 (u)) 

which from (P can be shown to satisfy 

[t(u),t(y)] = 0, \/u,veC. 

Above stro denotes the supertrace taken over the auxiliary space which is labelled by 0. 

In the usual manner the Hamiltonian associated with the transfer matrix is defined 
by the relation 



du 



M=0 
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Assuming regularity of the i?-matrix; i.e. 

R(0) = P 

where P is the (Z 2 -graded) permutation operator, yields 



N-l 



Hgiobai — /J Hi(i+i) + Hni (3) 
1=1 

where the local two site Hamiltonians are given by 

(4) 



H = ^-PR(u) 
du 



u=0 



2. Solutions arising from quantum superalgebras 

Let G denote a simple Lie (super) algebra of rank r with generators {ei, fi, hi\ r l=0 and let 
a\ be its simple roots. Here we adopt the convention that in the distinguished root basis 
ao labels the unique odd simple root. The quantum (super) algebra U q {G) can be defined 
with the structure of a (Z 2 -graded) quasi-triangular Hopf algebra ^?J. We will not give 
the full defining relations of U q {G) here (see e.g |3|, |]]) but mention that U q {G) has a 
coproduct structure given by 

k{hi)=I®hi + hi®I, A(a) = a <g> q- hl/2 + q ht/2 <g> a , a = e h /,. (5) 

Let Ti be an irreducible representation (irrep) of U q {G) afforded by the irreducible 
module V(A) where A denotes the highest weight. Assume that the irrep n is aflinizable, 
i.e. it can be extended to an irrep of the corresponding quantum affine (super) algebra 
U q (G). Consider an operator R(x) G End(V" ® V), where x G C is the multiplicative 
spectral parameter. It has been shown by Jimbo |JJ and also Zhang et. al. 11 for the 
^2-graded case that a solution to the linear equations 

R(x)A(a) = A(a)R(x) , Va G U q (G), 

R(x) (xvr(e^) ® n(q^/ 2 ) + n(q h */ 2 ) ® 7r(e*)) 

= (xTr(e^) ® 7i(q h * /2 ) + vr(g- /l * /2 ) <g> 7r(e*)) R(x) (6) 

satisfies the Yang-Baxter equation in the tensor product module V ® V <8> V: 

Ri2(x)R 13 (xy)R 23 (y) = R 23 (y)R 13 (xy)R 12 (x). (7) 

In the above, A = T • A, with T the twist map defined by 

T{a ® b) = (-l) [a][b] b ® a , Va, 6 G E^(G) 

and 7r(e$), 7r(h^) are matrices for operators associated with the highest root ^ such that 
the representation extends to a loop representation of the untwisted affine extension of 
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U q (G). The multiplicative spectral parameter x can be transformed into an additive 
spectral parameter u by x = exp(w). The structure of the Jimbo equations is such that 
the solutions obtained must necessarily have the difference property; i.e. 

R(u, v) = R{u — v). 

Explicit solutions of the Jimbo equations can be computed by tensor product graph 
methods as discussed in |28|, |29], |30|. 



As consequence of the Jimbo equations, for any Cartan element h the operator hi on 
the periodic lattice commutes with the global Hamiltonian ([3]) and so we may arbitrarily 
add such terms to the Hamiltonian and still yield a solvable model. In the next section, 
it will be described how the resultant model may be derived directly from an i?-matrix 
without the difference property. 



3. Twisting construction 

Let (A, A, R) denote a quasitriangular Hopf (super) algebra where A and R denote the 
co-product and .R-matrix respectively. Suppose that there exists an element F G A <g> A 
such that 

(A®I)(F)=F 13 F 23 , 
(/(g) A) (F) = F 13 F U , 

Fi 2 Fi 3 F 2 3 = F 23 Fi 3 F 12 . (8) 

Then (A, A F , R F ) is also a quasitriangular Hopf (super) algebra with co-product and 
i?-matrix respectively given by 

A i? = F 12 AF 1 - 2 1 , R F = F 21 RF 12 1 . (9) 

Throughout we refer to F as a twistor. 

The result stated above is a little more general than that originally proposed by 
Reshetikhin and is due to Engeldinger and Kempf |22] . In the original work pT| Reshetikhin 
imposed the additional constraint 

F 12 F 21 = 1(g) I 

and in the case that (A, A, R) is an affine quantum (super) algebra Reshetikhin gave the 
example that F can be chosen to be 

F = exp (hi (g) hj — hj (g) hi) <pij (10) 

i<j 

where {hi} is a basis for the Cartan subalgebra of the affine quantum (super) algebra and 
the <pij, i < j are arbitrary complex parameters. However following the construction of 
Engeldinger and Kempf it is possible to choose 



F = exp y %2(H i ®H j )<l) ij (11) 
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which obviously gives a twistor dependent on more free parameters. For our purposes 
either approach may be adopted but we will choose the latter for convenience. Note also 
that it is also possible (and for the construction below essential) to extend the Cartan 
subalgebra by an additional central extension (not the usual central charge) h which will 
act as a scalar multiple of the identity operator in any representation. 

It is worth observing here that the class of twistors (O) above also qualify as Hopf 
algebra preserving twistors of the type described by Andrews and Cornwell [j3TJj in their 
work on relating non-standard quantum algebras to standard ones. In notation as above, 
suppose there exists F G A <8> A which satisfies the following relations 

F12F23 = F23F12 
(A®/)F = F 23 F 13 
(I®A)F = F 12 F 13 . 

Then (A, A F , R F ) is also a quasi-triangular Hopf algebra with A F , R F given by eq. (|9]) 
above. 

For a given element h of the Cartan subalgebra we now choose a twistor in the par- 
ticular form 

F = exp(r]h ® h + rh <8> h) 
where 77, r are arbitrary scalars. This gives the twisted algebraic i?-matrix the form 

R F = exp(r]h <g> h + rh <g> h)Rexp(—r]h <8> h — rh (g) h) 

which still satisfies the MYB equation 

R12R13R23 — -^23-^13-^12- 

Suppose again that 7r extends to a loop representation of the affine quantum super- 
algebra. We let 7T eg) n £g) 7T act on the above equation but let the central element h take 
different values /3, 7, 5 in each of the tensor product spaces. Letting R F {u 1 [3 ) ^) denote 
the image of R F under n cg> ix where ho acts as (3 in the first space and 7 in the second, 
this then yields the matrix solution 

R 12 (u - v, (3, i)R 13 (u, (3, 5)R 23 (v, 7, 5) = R 23 (v, 7, 5)R 13 (u, (3, 5)R l2 (u - v, (3, 7) (12) 

^From this solution we may construct the transfer matrix 

t(u, (3, 5) = str (R 0L (u, {3, 6)....Ro 2 (u, (3, 5)R m (u, (3, 5)) (13) 

which as a result of QT3) forms a commuting family in two variables 

[t(u,{3,S),t(v, 7 ,5)) = 0. (14) 

It then follows that one may diagonalize t(u, (3, 5) independently of both u and (3. Setting 
exp h = M and R(u) as the image of R under it ® tt we get 

R F (u,P,5) = Ml 5 Mf R(u)M^ & M 2 ~ tP . 
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The above solution R F (u, ft, 5) can be made regular by setting 5 = 0, ft = u. The local 
two-site Hamiltonian operator associated with this solution is given by 

H F = ^-PR F ( u ,ft = u,5 = 0) 

du u=0 

= H + rih - rh 2 (15) 

with H given by (f|) above. It immediately follows that the global Hamiltonians are 
related by 

L 

H global = H global + (t) ~ t) ^kj. 

1=1 



5. The supersymmetric U model 



The quantum superalgebra U q (sl(2\l)) has simple generators {eo, fo, ho, e±, fi, hi} associ- 
ated with the Cartan matrix 

' 1 
-1 2 



A 



This algebra admits a non-trivial one parameter family of four-dimensional representations 
which we label n given by 



(e ) = VM e 2 + Vi a + !] e 



t(/o) 

7T(/l ) 

vr(ei) 

7T(/H) 



a{e\ + e\) + {a + l){el + e 



e 2 e 3- 



Above the indices of the elementary matrices e* carry the Z 2 -grading (1) 
(3) = 1 and we employ the notation 



(4) 



(16) 
0, (2) = 



q — q 
q-q- 



Associated with this representation there is a solution of the MYB equation with the 
difference property which is obtained by solving Jimbo's equations. The problem of 
obtaining this solution was considered in [|l(| Q (see also |33|). Adopting the prescription 
detailed in the previous section, the following solution of ( |I2"D is obtained with the choice 

h = (a + 1)1 — h — khi 



where A; is a free variable. 
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R(u,(3,S) = 

/Rll o 



E>12 









m 







Dl4 

H 14 




RP 





2i 



#23 






Rli 







pl4 
H 32 










d14 
-K 4 1 













d21 
^12 










n 14 



r?2i 

^21 









r>22 
^22 









K>23 








p24 
rt 24 













e>23 
"32 






n>23 

n 41 











R24 

rt 42 



Rll o 

R\\ 



Rll 






n>31 
"31 







Rll 









n.33 
"33 







r>34 
"34 





Rll 










p34 
"43 









V o o 









p41 

"23 













E>42 
_rt 2 4 






p41 
"32 










Ril 



where the non-zero entries are given by 

[u — a] [u — a — 1] 



nl2 _ 

rt 19 — 



exp[( V -r)((3-S)} 



b 12 



R 13 
R 

* r>12 
K 21 

* j^13 



exp[(r,-T)((l-k)P-5)] 
exp[-(?7 - r)5] 



[u + a] [u + a + 1] 

u] [u — a — 1] 



u + a] [u + a + 1] 
u] [u — a — 1] 



[u + a] [u + a + 1] 
«][«-!] 



ft 
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exp[((l-^)77-r)(/3-5)] g 
exp[(/c7? - r)(/3 - <5)]g' 
exp[-r(/3 - 5)}q 



[u + a] [u + a + 1] 



[w + a] [u + a + 1] 
[a] [if — a — 1] 



[w + a] [u + a + 1] 
2 „ L a][a + 1] 



-^23 



* pl4 
rl 32 



p21 



* D 22 



[u + a] [u + a + 1] 

«pt(-*r/J + <(* - 1), + r)S\ q ^ f'> + f>j. 

[u + a\ [u + a + 1] 

- exp[(* - l)rl, + (r - fe,)^ W^±^M 

[iz + a] [u + a + 1J 

«p[(i,-r)0> + (*-!)«], MI"" 



L 22 



exp[(l-fc)(7 7 -r)(/3-5)] 



[u + a] [u + a + 1] 
[u - a - 1] 



[« + a + 1] 



p41 







n42 
"42 










p43 
"43 








p44 
"44 
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*Rtl = exp[(v ~ r)(k(3 + (k - 1)5)] 
Rl\ = expp-l)(77-r)5] 



p44 
ti, 



[u + a][u + a + 1] 

M 



^ = - e xp[(^ + (A;-l)r)(/3-(J)]g 
*R 2 i = ex P [(k V + (k-l)r)(P-5)] 
Rll = e W [r(k - 1)((3 - 5)]q 



[u + a + 1] 

[a] [u — a — 1] 



[u + a] [u + a + 1] 

2g g2a+l q— 2a— 1 g2u+l _|_ q 2u-l 

^gU+a g— m— g— u— a— l^J 



*i?23 = _ exp [^ + ((l_A;)r-r7)5]g- u+1 / 2 
*i? 2 ? = _ eX p[(A;r ? -r)/3 + (l-A;)r5]g M+1 / 2 



[u + a+1] 

[a]^ 2 [a + l] l ' 2 [u] 
[u + a] [u + a + 1] 
[a^fa+l] 1 / 2 ^] 



[u + a] [u + a + 1] 

RH = exp[(^-r)(/3-M)] [ *~;-^ 

[u + a J [u + a + 1J 

= exp[(7 7 -r)((l-fc)/3-M)] 



u] 2 



= exp[^-r)(/3-5)] 
i?34 = ex p[-A;(77 - r)5] 



u + a] [u + a + 1] 

U — OL — 1] 



[u + a + 1] 

M 



i$ = -exp[(^-fcr)(/3-(J)]g 
= expKl-fc^-fcr)^-*)] 
= exp[-A;r(/3 - 6)]q 



[u + a + 1] 

[a] [u — a — 1] 



[u + a\[u + a + 1] 

2g -l _ g 2a+l _ q -2a-l + g -2«+l _ g -2u-l 
q— u— a^qU+a+1 q — u — ot — 1 ^ 

a + 1] 



*Rll = exp[(l-k)r]f3-(r]-kT)5]q- u - 1/2 
*Rll = exp[((l-k)r]-r)P + kr5]q 
< = exp[(77-r)/3] 



[u + a + 1] 

a] 1 /2[ Q; + i]V2[ w ] 



u + a] [u + a + 1] 
,-i/2 M^a+lFM 
u + a\[u + a + 1] 

u][u- 1] 



i?g = exp[(l-A;)(7 7 -r)/3] 
i?43 = exp[A;(r] — r 



[u + a] [u + a + 1] 

N 



[u + a + 1] 
[u] 



[u + a + 1] 



44 
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= -exp[ V (l-k)(/3-5)]q 
*< = -ex P {k V ((3-5)]q 



[u + a] [u + a + 1] 
« [«+!] 



[it + a + 1] 
[u + a + 1] 

- exp[(,-^ + (*-l), % ^W^±^ 

-, + i/2 M' /2 [" + i]' /2 N 

[it + a] [it + a + 1] 



*i$ = -exp[(ri + (k-l)T)P-kTi8\q- 



We remark again that the above i?-matrix solves the MYB equation subject to the 
rule (Q) which is a natural consequence of the superalgebra structure underlying the R- 
matrix. However, a change of signs in those matrix elements with a * gives a solution 
which satisfies the Yang-Baxter equation in the usual matrix form. 

The four dimensional space of states associated with the U q (sl (2\1)) representation 
may be identified with the electronic states 

|o>, IT), II), |U). 

Taking the logarithmic derivative of the transfer matrix yields the following global Hamil- 
tonian on a periodic lattice (with convenient normalization) 



H = - ^2(c\ a c {l+1)(7 + h.c.) exp[-l/2(r? - ae)n i{ _ a) - 1/2(77 + 

i.a 

+j ^( c !t c !i c (;+ih c (*+i)t + h - c -) + u ^2 n ^ n a 

+ <^ +t — - ( "- T)( tr~^' ,) )£'- 



n 



(*+i)(-ff)J 



(7; - r)(l - 2k)(q a+1 - q- a ~^ 



21ng 



where 



exp e = q 1 

la 



exp(-?7) = J^i, [/ = 2H" 1 



and the standard notation has been used for the Fermi and spin operators. The above 
model is the supersymmetric U model with chemical potetntial and applied magnetic field 
whose couplings may be chosen arbitrarily through the parameters 77 — r and k. 
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6. The Bariev model 



The Bariev model on a one- dimensional periodic lattice has the Hamiltonian 
H = - ^2(c] a c ii+1)a + h.c.) exp[-l/2(?7 - ar])n i{ - a) - 1/2(77 + 

It is apparent that the Bariev model is obtainable from the supersymmetric U model 
above in the limit a — > 00 by choosing 77 — r = 41ng and k — 1/2. It would thus appear 
reasonable to suspect that an i?-matrix solution for the Bariev model is obtainable from 
the solution for the supersymmetric U model. However, if one is to naively take the limit 
a — > 00 one finds 

lim R(u,P,5) = Ml {& - p) M'f- S) P 

a— »oo 

and as such the Bariev model is not obtainable in this fashion. 

Alternatively, we may make the following rescaling of the parameters 

u — > era, v — > era, 

r] — > 2i] In q, r — > 2r In g, 

and set 

/3 = u, 6 = v, k — 1/2. 
In this manner we obtain a solution of the MYB equation 

R 12 (u,v)R 13 (u,w)R 23 (v,w) = R 23 (v,w)R 13 (u,w)R 12 (u,v) 

in the limit a — > 00 where the matrix elements of t> ) now read 

oil _ 2(t,-t)(u-«) [ M ~ U ~ l ? 

n u — q 



Rll = q 



u — v + l] 2 

(r,-r)(u-2v) [U-V][U-V-1] 



[u-v + 1] 



2 



pl3 _ Jy-r)(u-2v) [u - v][u - V - 1] 

13 — y r i 112 

Dl4 _ -2fa-r)t; [U-V] 2 

n u — y i 



[u-v + l] 2 

*Dl2 _ ( t? -2t+1)(m-d) [ M ~ V ~ 1 ] 

21 ~~ y r i 112 

[ix — V + l\ z 

*Dl3 _ „(v-2t+1)(u-v) [u-V-l] 

3i ~~ y r i 112 

d14 _ (2-2r)(u-tQ 1 



*pl4 _ „(l-r)n+(-7?+2T-l)r) F V 
23 ~~ y T , -119 
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* pl4 
rt 32 




(1-t)u+(2t- V -1)v i U ~ V \ 

q [u - V + l] 2 


D21 
U 2l 




(v -r)(2u-v) [U-V][U-V-1] 

q [u-v + 1] 2 


* p22 
n 22 




(v-r)(u-v) [u — V — 1] 

q [u-v + 1] 


* p23 
-""23 




Jv-r)(u-v) [ U ~ V f 

q [u - V + l] 2 


p24 
n 24 




„-(v~t)v \ U - V \ 
CI 

H [u-v + 1] 


K 12 




(2n-T-l)(u-v) \U — V — 1\ 

q [u-v + 1] 2 


-^32 




n (V—T)(u—v) ' 

q [u-v + 1] 2 


p24 

n 42 




(1-t)(u-v) ' 

q [u-v + 1] 


* p23 

n u 




(t)-1)u+(t-2ti+1)v \- U ~ V ] 

q [u-v + 1] 2 


* z~>23 
K 41 




(r,-2T+l)u+(T-l)v \ u ~ v ) 

q [u-v + 1] 2 


p31 
^31 




(v -r)(2u-v) [U~V][U-V-1] 

q [u-v + 1] 2 


* d32 
-^32 




(n-T)(u-v) [u-v] 2 
q [u-v + 1] 2 


* d33 
-^33 




n (n-T)(u-v) [u-v-1] 
H [u-v + 1] 


D 34 
-^34 




n -(n-r)v W - v \ 
q [u-v + 1] 


d31 
K 13 




(2n-T-l)(u-v) [U — V — 1] 

q [u-v + 1] 2 


* d32 
-^23 




q [u-v + 1] 2 


d34 
-^43 




(1-t)(u-v) 1 

q [u-v + 1] 


* d32 
K 14 




„(v-l)u+T-2n+l)v l U ~~ V ) 

ff \ 1 1 1 V / 1 / 

q [u-v + 1] 2 


* d32 

-k 4 i 




(1+V-2t)u+(t-1)v i U ~ V \ 

q [u-v + 1] 2 


p41 




2( V -r)u W ~ V ? 

q [u-v + 1] 2 
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D42 _ ( V -r)u [ U ~ V \ 

rt 42 — y 



U — V + 1] 

43 = <? V 



^43 = („-r)t. [U-V] 



[u-v + 1] 



-R44 — 1 

d41 _ (2r? 

- rt i4 ~~ y 

-n, 24 — y 



* R 43 — 
-^34 - ~Q 



[u- 


■v + 


I] 2 


l)(u-v) 


1 




[u 


— v - 


f 1] 


l)(u-v) 


1 





-^23 - Q 



[u-v + 1] 

(2ri-T-i)u+(l-n)v l U ~ v \ 

u — V + l] 2 

[u — v] 



* p41 (2ry— t— r{)v 

-"■32 — y 



However, calculating the Hamiltonian associated with this solution one finds 
H = ~ I>Ui+l) CT + h - c -) + [V2(? - ^ X )(r - 77) + (g + y' 1 )] ^ 



which of course is simply a free fermion model with chemical potential and corresponds 
to a particular case of the Bariev model. 

7. Conclusions 

Here it has been demonstrated how i?-matrices which satisfy the MYB equation such that 
they do not have the difference property may be obtained. As an example, the anisotropic 
super symmetric U model with arbitrary chemical potential and magnetic field terms was 
derived, which, in a particular limit, reproduces the Bariev model. Regrettably, this 
approach does not shine a light on the connection with the i?-matrix solutions for the 
Bariev model computed in [T(| [□]] without the difference property. It is of course possible 
that a different limiting procedure may yield an i?-matrix solution for the Bariev model 
but its form is not apparent at present. 

An open problem for future work is to investigate whether Shastry's solution 0] and 
the recently introduced SU(N) Hubbard models [0, [H| may be related to an i?-matrix 



with the difference property via a twisting procedure. 
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